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Description 

[0001] The present invention pertains to the field of pulse monitors, and, more particularly, to the field of pulse mon- 
itors that use EKG signals to detect the pulse rate. 

5 [0002] In the prior art, two basic types of pulse rate monitors exist. The first type uses visible or infrared radiation 
which is projected through the skin to detect from radiation reflected from or penetrated through capillaries under the 
skin, pulsations of blood flow. Typically, these devices come in the form of a digital watch with a photodetector on the 
face of the watch or as a desktop unit or unit that clips to a belt with a clip connected to base unit, the clip for attaching 
to a fingertip or earlobe. Visible or infrared light passing through the skin is detected by the photodetector and gives 

10 an indication of pulsations in blood flow in capillaries. From these pulsations, the pulse rate is calculated. 

[0003] Numerous examples of this type pulse monitor exist as they are commonly found on fitness equipment such 
as treadmills, stationary bicycles and stair climbing machines. Another example of this type system is a pulse rate 
monitor wrist watch made by Casio. This watch reads both blood pressure and pulse. The watch has two sensors on 
top of the watch. The sensor on the lower left of the face of the watch is a photosensor which is to be covered with the 

is wearer's right index finger. The other sensor is to be covered by the right middle finger and is an electrode to pick up 
the EKG signals. The bottom plate of the watch body serves as the other electrode. 

[0004] There are several disadvantages to the photosensor/flow pulsation detectors. First, the finger position on the 
photosensor must be stable. Also, the force pressing the finger or ear lobe to the photodetector must be nominal. If 
the force is too high, the blood flow will be cut off, and no detection of blood flow pulsations can occur. If the force is 

20 too low, then any slight motion between the body and the sensor may cause inaccurate readings. Also, the reliability 
of the readings depends upon the ambient illumination (unless a separate photodiode supplies light for transmission 
through the skin) and upon the wavelength. Further, the flow pulse in a capillary looks like a sinusoidal waveform in 
shape. This makes it difficult to distinguish legitimate flow pulse signals from sinusoidally shaped noise waveforms. 
[0005] The second type of pulse rate monitor is the EKG type. These type pulse monitors work by picking up an EKG 

25 signal from the heart muscle itself and calculating the pulse rate from the EKG signal. One example of this type system 
is a pulse rate monitor wrist watch made by Casio. This watch reads both blood pressure and pulse. The watch has 
two sensors on top of the watch. The sensor on the lower left of the face of the watch is a photosensor which is to be 
covered with the wearer's right index finger and which detects fluctuations in light passing through the finger from 
ambient sources to determine when blood flow pulses occur. The sensor on the lower right of the face is one of two 

30 sensors for an EKG signal. The bottom surface of the watch is the other sensor for picking up the EKG signal. To use 
this watch to read blood pressure and pulse, the user must first use another independent instrument to measure blood 
pressure and pulse in an at-rest condition. These readings are then input to the watch. After inputting this data, the 
watch takes about 10-20 heartbeats with the user's fingers in contact with the two contacts on the face of the watch. 
During these 10-20 heartbeats, the watch learns the timing between the EKG signals picked up by the EKG contacts 

35 and the corresponding blood flow pulses in an at-rest condition. This timing serves as a reference for determination of 
blood pressure. The operative principle is determination by the watch of the timing between the EKG signal that causes 
the left ventricle to pump blood to the resulting pulse of blood flow detected by the photodetector in the capillaries of 
the wearer's fingertip. 

[0006] A drawback of this design is that the timing between the EKG pulse and the blood flow pulse changes with 

40 fitness level as the aerobic effect takes over and new biood flow paths are formed in the body. As a result, the manu- 
facturer recommends that the basic at-rest data read from an independent instrument be updated every three months. 
This is inconvenient unless the owner of the watch also owns independent Instruments to measure biood pressure and 
pulse rate. Further, the watch is incapable of measuring only pulse rate without aiso measuring blood pressure. 
[0007] EMG noise is particularly troublesome in EKG pulse monitors because its frequency is in the same range as 

45 the frequency of the sought after EKG signal. Therefore, special signal processing must be accomplished to separate 
EMG noise from the desired EKG signal. One type of signal processing methodology that has been tried in the prior 
art is autocorrelation. The EKG signal Is sensed by probes which are typically mounted on the handles of exercise 
equipment which the user grabs while exercising. The signals sensed by these probes, which contain the EKG signal, 
are passed through an autocorrelator which performs a correlation calculation between a piece of the signal represented 

50 by one buffer's worth of digital samples and an adjacent portion of the signal in time represented by another set of 
samples. Signal indication logic monitors the output of the autocorrelator for the presence of a periodic signal and 
generates a synthetic candidate heart rate signal that has the same frequency as the periodic signal in the output of 
the autocorrelator. The difficulty with this approach is that the EMG signals during exercise are also periodic and will 
cause peaks in the autocorrelator output, that do not represent periodic EKG signal. 

55 [0008] The class of products described above are generally used by health conscious people while working out. 
Because of the drawbacks of the approaches described above, there has arisen a need for a reliable pulse monitor 
that is small and can be worn on the wrist or attached to exercise equipment and requires no electrode gel or moisture 
or chest strap and which can accurately find the EKG signal despite low signal to noise ratio. 
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[0009] US-A-4938228 discloses apparatus according to the preambie of Claim 1 . 
[0010] The invention is set out in Claim 1 

A convenient low-cost heart rate monitor may be provided by embodiments of the present invention. In one embodiment, 
a digital filter structure Includes a low pass filter having a notch at 60 Hz and a bandpass filter which amplifies signals 
5 in a frequency range from 10-40 Hz and has a notch at 60 Hz, This digital filter has a recursive structure and uses 
integer coefficients to simplify and speed up the calculations. A four bit microcontroller may implement the digital filter. 
The output of the digital filter is subject to enhancement signal processing to emphasize QRS complexes indicative of 
human heartbeats. 

[0011] One embodiment of a heart rate monitor may include at least two electrical contacts for detecting the electrical 
10 signals when placed in contact with the body, an analog circuit that conditions the electrical signals, an analog-to-digital 
converter coupled to receive an analog output signal from the analog circuit and convert the analog signals to a plurality 
of digital samples, a digital filter for receiving the digital samples and suppressing noise signals that have frequencies 
below about 5-15 Hertz, and signals having frequencies above about 25-40 Hertz to generate filtered data, the digital 
filter being a recursive filter having integer coefficients, and an enhancement signal processor to receive the filtered 
15 data and highlight signals therein that have predetermined characteristics of QRS complexes in human heartbeat 
signals so as to generate enhanced digital data. 

[0012] In accordance with a method of an embodiment for distinguishing pulses caused by heartbeats from other 
pulses the method may include steps of discarding pulses that arrive at times inconsistent with plausible beat-to-beat 
variation in heartbeat rate as representing artifacts rather than heartbeats, and discarding pulses having amplitudes 

20 varying more than a predetermined amount from previously received pulses representing heartbeats. 

[0013] A further understanding of the nature and advantages of at least preferred embodiments of the invention 
herein may be realized by reference to the remaining portions of the specification and the attached drawings. 
[0014] Fig. 1 shows a typical application in a digital wrist watch for the EKG type heart rate monitor according to the 
teachings of the invention. 

25 [0015] Fig. 2 is a block diagram of the preferred embodiment of the electronics of the invention. 

[0016] Fig. 3 is an illustration of some of the waveforms involved in illustrating some of the issues dealt with by the 
teachings of the invention. 

[0017] Fig. 4 Is an overview flow chart of the software architecture according to the teachings of the Invention. 
[0018] Fig. 5 is a graph of the filter response of a digital low pass filter used in the preferred embodiment of the 
30 invention. 

[0019] Fig. 6 is a graph of the filter response of a digital bandpass filter used in the preferred embodiment of the 
invention. 

[0020] Fig. 7 Is a graph of the total filter response resulting from combining the effects of the digital low pass filter 
and the digital bandpass filter. 

35 [0021] Fig. 8 is a more detailed flow diagram of the enhancement signal processing of the right branch in Fig. 4. 
[0022] Fig. 9 is an illustration of the waveforms that illustrate the operation of the processing of Fig. 8. 
[0023] Fig. 10 is a flowchart describing steps of post-processing to determine the user's heart rate. 
[0024] Fig. 11 is a flowchart describing further steps of post-processing to determine the user's heart rate. 
[0025] Referring to Fig. 1 , there is shown a typical application, in a digital wrist watch for the EKG type heart rate 

*o monitor according to an embodiment The digital watch 1 0 includes a display such as LCD display 12 which is used 
to display the time, and when the pulse monitor mode is selected, displays the user's heart rate. In some embodiments, 
a separate display can be used to show the pulse rate with the time simultaneously displayed in display 12. The digital 
watch embodiment uses a three contact approach to help eliminate noise. Two electrical contacts 1 4 and 1 6 are placed 
on the front of the watch for the user to place his or her fingers on when the pulse mode is entered. A third electrical 

<s contact is located on the back side of the watch and is indicated in phantom at 1 8. The three contacts are connected 
to a differential amplifier inside the watch so that common mode noise is suppressed. It should be understood that the 
present invention is not limited to three contacts. The circuitry and software described below is incorporated into the 
watch 10 so that the unit is stand alone. 

[0026] Referring to Fig. 2, there is shown a block diagram of the preferred embodiment. Block 58 represents the 
so electrical contacts of any of the embodiments discussed above. Block 60 is a differential amplifier or instrumentation 
amplifier that is coupled via bus 62 to the contacts. The differential amplifier amplifies the analog signals on the three 
conductors of bus 62. The common Input on bus 62 is coupled to analog ground and the two remaining conductors 
are coupled to the plus and minus inputs of the differential amplifier. The differential amplifier serves to provide gain 
and to simultaneously eliminate common mode noise in the signal such 60/50 Hertz hum etc. The gain of the differential 
55 amplifier is set at a relatively low figure of 1-10 to prevent saturation of the operational amplifiers therein by the low 
frequencies that are still in the analog signal. 

[0027] Next, the output of the amplifier is filtered by an active analog bandpass filter 64. Active filters using RC 
components in operational amplifier circuits are used to implement the bandpass filter. The passband of this filter is 
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centered on some frequency between 5 and 20 Hertz, preferably 10-15 Hertz, and has a passband of approximately 
5-40 Hertz. The purpose of this bandpass filter is to remove high frequencies from the analog signal to prevent aliasing. 
The low frequency comer of the passband also eliminates DC components and any low frequency drift of the baseline 
caused by respiration muscle contractions or other muscle contractions such as typically occur when the user is ex- 

5 ercising. The gain of the bandpass filter is set between 1 and 10 at the center frequency. A passive bandpass filter 
could be used, but generally is not preferred since it requires too many components and the rolloff is not sharp enough. 
In the preferred embodiment, the passband filter is comprised of two separate hardware filters: a low pass filter with 
a first order rolloff and a comer frequency between 25-40 Hertz; and a high pass filter with a second order rolloff with 
a corner frequency of from 5-1 5 Hertz. The reason these two filters have different order rolloffs is to allow the differen- 

10 tiator to emphasize the high frequency components in the signal in a manner described below. In other words, the 
steeper rolloff on the low frequency corner provides better selectivity, while the less steep rolloff on the high frequency 
corner allows more high frequency component to reach the differentiator such that the differentiator will generate bigger 
slope numbers because of the sharp peaks in the signal represented by the data reaching the differentiator because 
of the higher content of high frequency components in the analog signals reaching the analog-to-digital converter. 

>5 These sharper comers cause the differentiation operation in the digital signal processing portions of the circuit to be 
described below to generate higher peaks. These higher peaks are filtered out by the moving average function, but in 
the process, the average rises thereby causing the EKG signals to stand out better against the background noise. 
[0028] Next, the analog signal is amplified in amplifier 66 which has a gain of 50-1000 so that the overall gain is 
about 1000-10,000. 

20 [0029] Next, the analog output of the amplifier 66 is applied to the input of an analog-to-digital converter (hereafter 
AD converter) 68, which, in the preferred embodiment is integrated onto the microcontroller integrated circuit. The 
analog signal is converted to 8-bit digital samples at a sampling rate of 1 80 Hertz, or some other multiple of 60 Hertz 
preferably. Other sampling rates can also be used so long as the high frequency components of the real signal from 
the electrodes have been filtered out prior to the real signal reaching the analog-to-digital converter. If these high 

25 frequency components have not been filtered out, alias signals will appear in the EKG frequency range resulting from 
beating of the sampling rate with the high frequency components. This is why use of an analog bandpass filter in front 
of the analog-to-digital converter is preferred. Different numbers of bits in each sample could be used as well as a 
higher or lower sampling rate, but 1 80 Hertz is preferred to make the digital filtering easier since the 1 80 Hertz sampling 
rate is a multiple of the 60 Hertz frequency of a common source of noise. Further, the coefficients for the digital signal 

30 processing stages described below are set for a sample rate of 1 80 samples per second, so If other sample rates are 
used, the coefficients must be optimized again. In some embodiments, the AD converter is part of a microcontroller 
integrated circuit 70. The AD converter should be able to operate on a 3 volt supply or some supply voltage easily 
achievable with batteries and, for size considerations In some embodiments, preferably has a serial format output data 
stream if the AD converter Is not integrated onboard the microcontroller. 

35 [0030] The digital data stream from the AD converter is input to the microcontroller 70 for further signal processing 
to be described below. The microcontroller is preferably a 4-bit or 8-bit machine such as a Samsung KS57C2408 or 
Hitachi H8/3612 or equivalent. These microcontrollers have on-board analog-to-digital converter, dual clocks, LCD 
drive and buzzer drives. There Is no need for an on-board DSP or a separate DSP since the digital filter coefficients 
are integers which substantially simplifies and speeds up processing in the software filtering modules. The microcon- 

<o troller is also coupled to a number of user interface peripherals represented by block 72. The user interface peripherals 
include LCD displays, mode control switches, buzzer(s), sensors for monitoring track/wheel speed or as part of an 
automatic feedback loop which will cause more drag on the wheel/track when the heart rate drops below a target range, 
and output ports to which different resistors may be connected to change the gain of the operational amplifiers, output 
signals to implement automatic load variation in embodiments where the drag on the belt/wheel is automatically varied 

45 in accordance with the user's heart rate, such as circuitry to change the gears on a mountain bike to a lower gear or 
higher gear automatically based upon the user's heart rate, and alarm units to give users audible alarms when some- 
thing is not working or needs adjustment. The microcontroller Incorporates a main clock 74 which can be shut down 
in power down mode while an independent digital clock also incorporated within the microcontroller is not disabled 
during power down mode. 

so [0031] In some embodiments, a separate digital clock integrated circuit 76 with its own clock 78 serves to keep track 
of the time/date etc. and control the display thereof. In alternative embodiments, the digital clock chip can be any digital 
clock chip where display 80 displays only time and date so long as the clock chip has an interface to the microcontroller 
70 to receive heart rate data. However, in the preferred embodiment, digital clock IC 76 is Incorporated Into the micro- 
controller and receives for display data encoding the heart rate reading, calories burned, speed, distance covered, 

ss exercise time etc. via the conductors of bus 82 for display on LCD display 80. 

[0032] There follows a discussion of the digital signal processing implemented by the firmware which controls oper- 
ations of the microcontroller 70. However, to better understand the approach used, the characteristics of a typical EKG 
signal and an EKG signal with noise will first be discussed. Referring to Fig. 3, there is shown a number of waveforms 
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on time lines (A) through (G). Time line (A) shows a typical EKG signal, which is sometimes referred to as a QRS 
complex. An EKG signal is comprised of a first small positive peak designated as P on time line (A), followed by a 
negative peak designated Q. After the negative Q peak, a very large positive peak immediately follows which is des- 
ignated as R on time line (A). The R peak Is Immediately followed by a negative S peak which is immediately followed 

s by a smaller positive T peak. The slope of the R peak is symmetrical on both sides and has a unique characteristic value. 
[0033] Time line (B) shows a EKG signal superimposed with 60 Hertz line interference. Time line (C) shows the result 
of passing the signal of time line (B) through a second order Butterworth lowpass filter. The signal on time line (C) still 
exhibits some 60 Hertz noise. Time line (E) represents and EKG signal superimposed upon a low frequency noise 
sources such as respiration EMG or exercise EMG signals. Time line (F) shows the result of passing the signal of time 

io line (B) through a second order Butterworth high pass filter showing distortion of the signal from the action of the filter. 
Time line (D) shows the signal resulting from passing the noisy signal of time line (B) through an integer coefficient low 
pass filter such as is used in the invention. Time line (G) shows the signal resulting from passing the signal of time line 
(E) through an integer coefficient high pass filter such as is used in the invention. Both low and high pass integer 
coefficient filters are used in the invention to eliminate both 60 Hertz noise and low frequency noise. 

15 [0034] Referring to Fig. 4, there is shown a flow diagram of the main processing carried out by microcontroller 70 
on the digital samples received from the AD converter 68. The digital filtering which results in the signals on time lines 
(D) and (G) in Fig. 3 is represented by the digital filtering block 90. The function of this block is to remove high and low 
frequency noise from the digital samples. 

[0035] The digital filtering steps represented by block 90 comprises, in the preferred embodiment, a first step of low 
20 pass filtering the incoming digital data with a low pass filter having a notch at 60 Hertz and then passing the resulting 
filtered data through a bandpass filter which amplifies the signals in a frequency range from 10-40 Hertz and has a 
notch at 60 Hertz. The reason a low pass filtering step is used prior to passing the data through is to remove remaining 
undesired high frequency components in the incoming data prior to passing the signals through the bandpass filter for 
amplification and further filtering. 
25 [0036] Fig. 5 represents the frequency response of the preferred embodiment for the digital low pass filter imple- 
mented by the preferred embodiment of the digital filtering step 90 in Fig. 4. This digital filtering is achieved by imple- 
menting a recursive filter having integer coefficients to simplify and speed up the calculations required to implement 
the filter. However, any digital or analog filter that implements a frequency response roughly equivalent to the frequency 
response of Fig. 5 will suffice for purposes of practicing the invention, e.g., non-recursive, recursive with different integer 
30 coefficients or floating point coefficients, separate DSP etc., analog or any combination of the above. In fact, the low 
pass filtering step can be eliminated altogether in some embodiments. The purpose of the low pass digital filter is to 
remove 60 Hertz noise and to enhance signal to noise ratio. 

[0037] Fig. 6 represents the frequency response of the preferred embodiment of the bandpass filter implemented as 
the second stage of digital filtering represented by block 90 in the preferred embodiment. This digital filtering frequency 

55 response is achieved by implementing a recursive filter having integer coefficients to simplify and speed up the calcu- 
lations required to implement the filter. However, any digital filter that implements a frequency response roughly equiv- 
alent to the frequency response of Fig. 6 will suffice for purposes of practicing the invention, e.g., non-recursive, re- 
cursive with different integer coefficients or floating point coefficients, separate DSP etc. or any combination of the 
above. The bandpass filter function could also be done using an analog filter One purpose of the bandpass filter is to 

40 reduce low and high frequency EMG signals caused by exercise, respiration and muscle tremor in addition to further 
reducing 60 Hertz hum noise. Some of the above mentioned EMG signals lie outside the passband and are therefore 
attenuated. EMG and other noise signals In the passband are removed by a learning, adaptive threshold detect oper- 
ation and further uses rules of reason implemented in the software to arbitrate whether pulses which exceed the adap- 
tive threshold are true heart beats or are artifacts or noise. 

45 [0038] The low pass filtering step can be eliminated in some alternative embodiments since the bandpass filter also 
removes the high frequencies. However, the low pass filter is preferred because if it is not present, there is less effective 
filtering out of 50 Hertz noise which is only an Issue If the invention is to be used in any European or other country 
where 50 Hertz power is used. 

[0039] Retumingto the consideration of the overall digital signal processing represented by Fig. 4, after digital filtering 
so is performed, two different processes are performed to find the EKG signal in the filtered data. In the preferred embod- 
iment, only the enhancement signal processing represented by block 92 is performed. In another embodiment, template 
matching or cross-correlation Is also used to provide alternative results which can be cross-checked against the results 
of the enhancement signal processing. 

[0040] Referring to Fig. 8, there is shown a flow chart of the basic operations that are performed in the enhancement 
55 signal processing block 92. The first step is differentiation, represented by block 96. The reason this differentiation 
operation is performed is because the filtering process removed some of the information content of the EKG signal, 
but there Is still a fairly high spectral content of high frequency Fourier components of the EKG signal caused by the 
Q, R and S peaks which is within the passband of both the analog and the digital filters. To enhance this component, 
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the differentiation process is performed on the array of sample data to find the slopes of all peaks in the data, The 
incoming filtered data has a signal waveform which can look like the waveform at time line A in Fig. 9. The differentiation 
process receives this waveform and outputs a waveform that looks like the signal on time line B in Fig. 9. 
[0041] The slopes of the two sides of each peak in the signal on time line A are successively positive and then 

5 negative. Because a moving average (an integration step) is to be calculated to smooth out the output of the enhance- 
ment signal processing output and to improve the signal to noise ratio, negative peaks are undesirable because they 
tend to cancel positive peaks and make the signal to noise ratio worse. To get rid of the negative peaks and possibly 
improve the signal to noise ratio, in the preferred embodiment, a squaring operation represented by block 98 is per- 
formed. The squaring operation outputs a waveform like that shown on time line C of Fig. 9. Note that ail negative 

10 peaks are converted to positive peaks. 

[0042] Because squaring requires multiplication, and because multiplication is slow and consumes too much com- 
puting resource, the preferred embodiment is to substitute a look up table containing the squares of all values likely to 
be found in the input signal to avoid the processing necessary in multiplication. This is much faster than actually per- 
forming the multiplication. In an alternative embodiment, an "absolute value" step can be performed for block 98. This 

15 absolute value operation receives both positive and negative peaks and outputs a signal which is all positive peaks, 
each representing the absolute value of the magnitude of a corresponding positive or negative peak. Squaring is pre- 
ferred since this has the effect of multiplying the size of the peaks in the signal thereby enhancing the signal to noise 
ratio. More than one squaring operation may be applied, preferably two, to further enhance the signal-to-noise ratio. 
[0043] Finally, a moving average is computed as represented by block 100. The purpose of calculating this moving 

20 average is to smooth out any sharp spike-like peaks in the output of the squaring operation or absolute value step. 
The moving average computation is like an integration which calculates the area under each peak in the output from 
the squaring or absolute value step. The moving average step also increases the signal to noise ratio by taking each 
pair of peaks in the signal output by block 98 and summing the area under each curve. This enhances the amplitude 
of the resulting single peak which results from each pair of adjacent peaks in the signal on time line C. For example, 

25 peak 102 results in the output of the moving average step from the integration or filtering of the pair of peaks 104 and 
106 that arrived from the squaring step. The moving average is calculated by taking a predetermined number of se- 
quential samples and adding their values and then dividing the sum by the number of samples so added. The individual 
samples in each moving average can be weighted in some embodiments. This process is carried out continuously for 
each contiguous group of samples. 

so [0044] Equation (1) below represents the mathematical expression which completely specifies the digital signal 
processing which implements the recursive low pass filter step of block 90 in Fig.- 4 with a notch at 60 Hertz. 



Equation (1): Y n = 1/8*(2Y n . 1 - Y n . 2 + X n - 2 + X^) - 

35 

The Y terms represent previous outputs such that Y rv1 represents the most recent previous output while Y n . 2 represents 
the second most recent previous output. The X terms represent recent data, except that the result of a previous stage's 
calculation, e.g., Y n in equation 1 becomes raw data point X n for the next stage, e.g., X n in Equation (2) below is actually 
the result Y n from the calculation process represented by Equation (1). For example X n represents the most recent 

40 data while X^ represent the fourth most recent data and so on. Note how the coefficients of the recursive filter spec- 
ification represented by Equation (1 ) are all integers thereby greatly speeding up and simplifying the calculation. There 
are many different combination of coefficients which can result in acceptable filter transfer functions. There may even 
be other filter specifications with more points or different subscripts which result in acceptable filter responses. Equation 
(1) only represents one acceptable example, and any alternative that results in a similar range of frequencies in the 

45 passband as shown In Fig. 5, and relatively similar gain levels will be acceptable. The factor 1/8 in the equation of 
Equation (1) is a normalization factor based upon the gain of the filter. 

[0045] Equation (2) represents the mathematical expression which completely specifies the digital signal processing 
which implements the recursive bandpass filter step of block 90 in Fig. 4 with a notch at 60 Hertz. 

50 Equation (2): Y n = 1/12*(2Y fVl - SY^ + 2Y n . 3 + X n - 2 X^ + X,, 12 ) 

The same can be the about other possible alternatives to the filter specification of Equation (2) as was said about 
Equation (1). The factor 1/12 in the equation of Equation (2) is a normalization factor based upon the gain of the filter. 
55 [0046] Equation (3) represents the mathematical expression which completely specifies the digital signal processing 
which Implements the differentiation step 96 in Fig. 8. 
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Equation (3): Y n = 1/4'(2' X n + X fr1 - X n _ 3 -2* X^) 

Many other combinations of weighting coefficients, numbers of sample points and possibly even different subscripts 
5 exist. For example, the denominator 4 in Equation (3) is chosen so as to make the result of the differentiation process 
larger such that the input data to the squaring process is larger. This results in greater dynamic range. However, other 
denominator factors exist which can be used without departing from the invention defined in the claims. Any such 
combination that result in a filter characteristic that has the same relative range of frequencies within the first major 
loop of the transfer function where most of the useful work gets done will be acceptable to practice the invention. For 
10 example, different numbers of sample points with a different divisor will also work to practice the invention, 

[0047] Equation (4) represents the mathematical expression which completely specifies the digital signal processing 
which implements the moving average step 100 in Fig. 8. 



15 Equation (4): 

20 

The moving average function has the effect of amplifying the low frequency components. The number of points selected 
for the -window" over which the moving average is calculated is important in that it defines where the first notch is in 
the transfer function which has the characteristic of defining how high the amplitude is of the resulting moving average 
peak. Eight points were selected in the preferred embodiment to enhance the signal-to-noise ratio by raising the am- 

25 plitude of the peak that results from each peak pair in the signal from the squaring or absolute value operation of block 
98 in Fig. 8. Other values for signal-to-noise ratio would also work, but perhaps not as well. For example, 16 points 
would smooth the output more, but the amplitude of the peak is lower. A number of points fewer than 8 may result in 
a sharper peak with more adjacent noise thereby resulting In lower signal to noise ratio. The 1/64 factor In front of the 
summation is a scaling factor which can be used as a digital gain control factor, In some embodiments, the adaptive 

30 learning process will be used to affect this digital gain control factor based upon characteristics found in the sample data. 
[0048] Referring again to Fig. 4, in an alternative embodiment, a second type of signal processing of the digital data 
after digital filtering is carried out as a cross check against the results of the enhancement signal processing or as an 
alternative thereto. The alternative or additional processing Is symbolized by block 1 08 and involves template matching 
in the preferred embodiment or cross-correlation where adequate processing power is available to do the multiplications 

35 involved in the cross-correlation calculation. The template matched filtering is basically a crude but fast way of per- 
forming rudimentary cross-correlation to determine the degree of likeness between two waveforms. Template matching 
is discussed in greater detail in the parent application, U.S. Patent App. No. 08/654,373. 

[0049] The processes of block 108 and 92 in Fig. 4 have the effect of removing artifacts and noise and enhancing 
the peaks in the incoming data which are likely to be actual QRS EKG waveforms. In the process of block 92, the 

40 differentiation step of block 96 in Fig. 8 tends to remove or suppress EMG peaks and enhance EKG peaks because 
the EKG peaks have steeper slopes and sharper tips than typical EMG waveforms. This results in larger numbers 
being output from the differentiator for EKG signals that have been differentiated than for EMG signals. This fact trans- 
lates to higher peaks in both the positive and negative amplitudes in the signals on time line B of Fig. 9 representing 
the differentiator results from differentiating EKG QRS waveforms than result from differentiation of typical EMG pulses. 

45 This further suppresses EMG noise and enhances the signal to noise ratio. 

[0050] The template matching or cross-correlation process of block 1 08 removes artifacts by virtue of suppressing 
peaks that do not have the shape of the template. The template is designed to simulate closely the shape of the QRS 
waveform in a typical EKG waveform. In some embodiments, the shape of the template will be adjusted "on-the-ffy" 
after actual EKG signals are isolated for this particular user so as to more closely match the shape of the QRS waveform 

so of this particular user. Suppressing of EMG waveforms and other noise and artifacts occurs because these spurious 
signals do not have the same shape as the template and therefore result in lower correlation peaks in the output data 
of the template matching process. 

[0051] After enhancement of the digital data and removal of as many artifacts as possible in the enhancement proc- 
ess, either by the process of block 92 or block 108 in Fig. 4, the resulting digital data is analyzed to determine the 
55 user's current heart rate. This is done by post-processing heart rate determination module 109 in Fig. 4. The process 
represented by block 1 09 could be any one of a large number of different types of discrimination or arbitration routines 
to decide which of the pulses in the filtered, enhanced data from either path 11 0 or path 112 in Fig. 4 represent actual 
heart beats and to calculate the heart rate therefrom. 
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[0052] Fig. 10 is a flowchart describing steps of post-processing to determine the user's heart rate in accordance 
with one embodiment of the present invention. Fig. 1 0 shows steps of first establishing that a series of pulses correspond 
to the user's heartbeat The steps of Fig. 10 rely on both the amplitudes and positions of pulses to establish that a 
series of pulses correspond to heartbeats rather than noise artifacts that have survived the previous processing stages. 

5 At step 1 002, post-processing heart rate determination module 1 09 awaits the first pulse above a predetermined am- 
plitude threshold. In one embodiment, this threshold is 1/4 of the full scale value of the A/D converter output. Once this 
pulse arrives, post-processing heart rate determination module 1 09, at step 1 004, looks for a second pulse having an 
amplitude within 50% of the amplitude of the first pulse and arriving at a time indicative of a realistic human heart rate, 
e.g., between 40 and 240 bpm. 

10 [0053] If no second pulse meeting these criteria arrives, processing returns to step 1002 to reawait a first pulse. It 
should be noted that if a second pulse arrives but does not meet the criteria for a second pulse, it may meet the criteria 
for a first pulse and step 1 002 will be skipped. If a second pulse meeting these criteria arrives, this suggests that both 
the first and second pulse may both be heartbeats, At step 1 006, post-processing heart rate determination module 1 09 
awaits a third pulse having 50% of the amplitude of the average amplitude of the first and second pulses and arriving 

'5 at a time defining an inteipulse period having a duration within 6.75% of the interval between the first putse and the 
second pulse. 

[0054] If no such third pulse arrives, the first two pulses are discarded and processing returns to step 1 002 to await 
a new first pulse. If a third pulse meeting the criteria arrives, post-processing heart rate determination module 109 
begins a process of confirming the current measured heart rate by awaiting a fourth pulse. However, post-processing 

20 heart rate determination module 1 09 takes extra precautions to insure that the fourth pulse is not noise. First, at step 
1 008, post-processing heart rate determination module 1 09 waits for 50% of an average period determined to be the 
average of the period between the first and the second pulses and the period between the second and the third pulses. 
[0055] Post-processing heart rate determination module 1 09 then looks for a noise pulse that might begin before the 
expected arrival of the fourth pulse and extend into the expected arrival period. In one embodiment, the threshold for 

25 detecting such noise pulses starts at approximately 50% of the mean amplitude of the first three pulses and decreases 
in four steps to approximately 19% of the mean amplitude. At step 1010, post-processing heart rate determination 
module 109 spends the next 1 2.5% (or fourth of the remaining time) of the period until the expected arrival of the fourth 
pulse checking for a noise pulse having an amplitude above the initial noise pulse threshold of 50%. If such a pulse 
begins but does not end during this subperiod, then at step 1012, post-processing heart rate determination module 

30 1 09 sets a noise pulse flag. If no such pulse arrives, or after step 1012, post-processing heart rate determination module 
1 09 checks if a previous noise pulse ends during this subperiod without a new one starting, at step 1 01 4. If a previous 
pulse ends, then the noise pulse flag clears at step 1016. After the noise pulse flag clears, or if the noise pulse flag is 
not cleared, at step 1018, the noise pulse threshold is lowered for the next subperiod of the last 50% of the' period 
before the expected arrival of the fourth pulse. Then post-processing heart rate determination module 109 checks at 

35 step 1020 whether, the next 12.5% subperiod is the last before the expected arrival of the fourth pulse. If the next 
subperiod is not be the last, post-processing heart rate determination module 109 returns to step 1010 to adjust the 
noise pulse flag. 

[0056] If the next 12.5% subperiod is the last, at step 1022, the preferred embodiment sets and clears the noise 
pulse flag based on any detected pulse activity. The noise pulse threshold has now stepped down to its minimum. 

40 Once the last 675% of the expected inteipulse period arrives, a new pulse may be the expected heartbeat. If, at step 
1 026, a pulse arrives within 6.75% of the expected arrival time, and having amplitude within 50% of the average of the 
first three qualified pulses, and while the noise pulse flag is clear, then this is considered to be the fourth qualified pulse 
and the heart rate based on the first four qualified pulses is read out at step 1 028. Requiring the noise pulse flag to be 
clear eliminates artifacts caused by earlier noise pulses extending into the expected arrival time of the fourth pulse. If 

4* no such fourth pulse arrives, the first three pulses are discarded and post-processing heart rate determination module 
109 proceeds to step 1002 to await a new first pulse. 

[0057] After step 1028, post-processing heart rate determination module 109 may be considered to be in a steady 
state, with a current displayed heart rate being based on the current average period between pulses. However, post- 
processing heart rate determination module 109 continues to guard against erroneous readings caused by noise. The 

so process of qualifying successive pulses as heartbeats is similar to the steps of qualifying the fourth pulse as shown in 
Fig. 10. In facts, steps 1102, 1104, 1106, 1108, 1110, 1112, and 1114 are essentially identical to steps 1002-1020 of 
Fig. 10. In Fig. 11 , post-processing heart rate determination module 109 looks for noise pulses during the last 50% of 
the period between the last heartbeat and the expected arrival time of the next heartbeat as determined by the mean 
interpulse period. At step 1116, during the last 12.5% of this interpulse period, post-processing heart rate determination 

55 module 109 continues to set and clear the noise pulse flag based on received noise pulses but also awaits the next 
heartbeat pulse during the last 6.75% of the interpulse period. At step 1118, post-processing heart rate determination 
module 1 09 determines if a pulse arrives within 50% of the median amplitude, within 6.75% of the average period, and 
while the noise pulse flag is clear. If such a pulse arrives, heartbeat display continues at step 1120. 
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[0058] If no such pulse arrives, then post-processing heart rate determination module 1 09 checks at step 1 1 22 wheth- 
er a qualified pulse has in fact arrived in the last five seconds. If no qualified pulse has arrived in the last five seconds, 
it is necessary to reacquire the heartbeat at step 1124 by repeating the steps of Fig. 10. If a qualified pulse has arrived 
In the last five seconds, heartbeat display continues at step 1120 and post-processing heart rate determination module 
109 repeats the steps of Fig. 11 , looking for a new qualified pulse at the next expected arrival time. 
[0059] In the foregoing specification, the invention has been described with reference to specific exemplary embod- 
iments thereof. It will, however, be evident that various modifications and changes may be made. 



Claims 

1 . Apparatus for determining heart rate from electrical signals generated within a body, comprising: 

at least two electrical contacts (58) for detecting said electrical signals when placed in contact with said body 
an analog circuit (60), (64), (66) that conditions said electrical signals; 

an analog-to-digital converter (68) coupled to receive an analog output signal from said analog circuit and 
convert said analog signals to a plurality of digital samples; 

a digital filter for receiving said digital samples and suppressing unwanted signals (70); 
an enhancement signal (92) processor to receive said filtered data and highlight signais therein that have 
predetermined characteristics or QRS complexes in human heartbeat signals so as to generate enhanced 
digital data, 

characterised In that said enhancement signal processor comprises: 

a differentiator (96) for determining the slope of peaks in said filtered data and generating a slope signal which 
defines the magnitude and sign of the slopes of each portion of each said peak; 
a squaring processor (98) for squaring the results from said differentiator; and 

a moving average processor (100) for computing a moving average of said positive values only signal and 
outputting a moving average signal which defines said moving average over time. 

2. The apparatus of claim 1 wherein said analog circuit comprises: 

a differential amplifier (60) having a gain of from approximately 1 -1 0 coupled to said electrical contact means 
for amplifying any signals detected by said electrical contact means and suppressing any common mode noise; 
and 

an analog bandpass filter (64) coupled to receive the output of said differential amplifier and having a passband 
from approximately 5 to 40 Hertz, said bandpass filter comprising a low pass analog active filter having a first 
order rolloff characteristic with a comer frequency between 25 and 40 Hertz, and a high pass analog active 
filter having a second order rolloff characteristic with a corner frequency between 5 and 15 Hertz. 

3. The apparatus of claim 1 or claim 2 wherein said analog-to-digital converter (68) uses a sample rate that is a 
multiple of the powerline AC voltage frequency; 

4. The apparatus of any preceding claim wherein said digital filter (90) receives said digital samples and suppressing 
noise signals that have frequencies below 5-15 Hertz and signals having frequencies above about 25-40 Hertz to 
generate filtered data, said digital filter being a recursive filter having Integer co-efficients. 

5. The apparatus of any preceding claim wherein said differentiator is a digital signal processor which computes the 
mathematical expression 

Y n = [2X n + X M .X^ 3 - 2X^1/4, 

where 

Y n = represents slope at any particular sample time N, 

2X n s represents twice the amplitude of the most recent digital data sample from the output of the previous filter, 
X n ..j m represents the amplitude of the next most recent digital data sample from the output of the previous filter, 
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X n . 3 = represents the amplitude of the third most recent digital data sample in the stream from the output of 
the previous filter, and 

2X n . 4 = represents twice the amplitude of the fourth most recent digital data sample from the output of the 
previous filter. 

5 

6. The apparatus of any preceding claim wherein said moving average processor calculates said moving average 
as a stream of digital data samples with the most recent digital data sample representing the moving average 
represented by Y N and wherein each Y N in said moving average data stream is calculated by summing the 8 most 
recent data samples in said stream of data samples representing said filtered data and dividing the sum by 64. 

10 

7. Apparatus according to any preceding claim, further comprising a post processing heart rate determination signal 
processor arranged to analyze said digital data to determine the heart rate therefrom. 

8. The apparatus of Claim 7 further comprising a display for displaying said heart rate determined by the post-process- 
15 ing heart rate determination signal processor means (1 09). 

9. Apparatus according to any of Claims 1 to 6, further comprising post-processing heartbeat determining means 
configured to distinguish heartbeat pulses from other pulses by virtue of being configured to: 

20 discard pulses that arrive at times inconsistent with plausible beat-to beat variation in heartbeat rate as rep- 

resenting artifacts rather than heartbeats; and 

discard pulses having amplitudes varying more than a predetermined amount from previously received pulses 
representing heartbeats. 

25 10. Apparatus according to Claim 9 further configured to; 

identify an artifact pulse arriving before an expected arrival time of a heartbeat pulse; and 

Identify a pulse that arrives at said expected arrival time as a heartbeat pulse only if said artifact pulse has 

subsided. 

30 

11. Apparatus according to Claim 9 or 10 configured to discard pulses that arrive at a time indicative of a change of 
period between heartbeats greater than approximately 6.75%. 

12. Apparatus according to Claim 9, 10 or 11 configured to discard pulses that have amplitudes varying more than 
35 50% from previously received pulses representing heartbeats. 

13. Apparatus according to Claim 10 configured to: 

after approximately 50% of a period between a last heartbeat pulse and said expected arrival time of a heart- 
<o beat pulse, set a noise pulse threshold for detecting artifact pulses; 

lower said noise pulse threshold until said expected arrival time; and 

Identify non-heartbeat pulses above said noise pulse threshold as artifact pulses. 

14. A method of detecting a heartbeat rate among unwanted signals, comprising the steps of: 

45 

sensing signals from a body using electrodes in contact with the skin of the body; 
amplifying the signals; 

filtering out noise below and above the frequency range in which the desired heartbeat rate will lie so as to 
provide a filtered analog signal; 
so converting the filtered analog signal to a plurality of digital samples; 

digitally filtering the digital samples to remove unwanted signals so as to generate filtered digital samples; 
digitally processing the filtered digital samples to enhance heartbeat peaks represented therein to provide 
enhanced digital samples, characterised by 

differentiating the filtered digital samples to determine the slope of peaks encoded therein to provide differen- 
55 tiatior digital samples, 

squaring the amplitude of signals encoded in said differentiator digital sample to provide squared digital sam- 
ples, and 

computing a moving average of said squared digital samples; and processing said enhanced digital samples 
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to determine said heartbeat rate. 

15. The method of Claim 14 wherein the step of digitally filtering the digital samples comprises the step of recursively 
filtering the sample data derived from the signals sensed from said body according to the following recursive filter 

5 specifications, recursively low pass filtering according to the equation: 

Y n = 1/8(2Y rH -Y n . 2 + X n -2X n . 4 + 

10 where, 

Y n represents the output digital sample on each iteration, 

Y fr1 and Y^ represent previous outputs such that Y rv1 represents the most recent previous output while Y n . 2 
represents the second most recent previous output, 

15 

and wherein 

X terms such as X n , X^ and X n . 8 represent recent raw data inputs such that X n represents the most recent 
raw data input, X n . 4 represents the 4th most recent raw data input and X^ represents the 8th most recent raw 
data input, and 

20 recursively bandpass filtering according to the equation: 

Y n = 1/1 2 (2Y^ - 9f^ + 2Y rv3 + X n - 2X..6 + X^) 

25 where, 

Y n represents the output digital sample of the bandpass filtering step on each iteration, and 
Y n .., and Y^ represent previous output digital samples of the bandpass filtering step such that Y n-1 represents 
the most recent previous output while Y n . 2 represents the second most recent previous output, and Y^ rep- 
30 resents the third most recent previous output, 

and wherein 

X terms such as X 0t X n _ 4 and X^q represent recent sample inputs such as samples derived from the raw or 
analog filtered data or digital sample results of a previous stage's calculation such as the digital low pass filtering 
35 step. 

16. The method of Claim 14 wherein said step of digital filtering removes further remnants of frequencies above and 
below the range of frequencies in which said heart beat will lie and to suppress powerline hum at approximately 
50-60 Hertz. 

40 

17. The method of Claim 16 wherein said step of computing said moving average comprises the steps of: 

using a computer to calculate a moving average of said squared digital samples according to the mathematical 
expression, 

45 




50 

where 

Y N represents moving average at any particular, sample time N, 
55 x N .| represents the amplitude of the digital data sample at time N-i. 

1 8. The method of Claim 1 4 wherein said squaring of the amplitude is performed by looking up squares in a lookup table. 



EP 0 861 045 B1 



1 9. A method according to any of Claims 1 3 to 1 8, further comprising distinguishing pulses caused by heartbeats from 
other pulses by the steps of: 

(a) discarding pulses that arrive at times inconsistent with plausible beat-to beat variation heartbeat rate as 
representing artifacts rather than heartbeats; and 

(b) discarding pulses having amplitudes varying more than a predetermined.amount from previously received 
pulses representing heartbeats. 

20. The method of Claim 19 further comprising the steps of: 

(c) identifying an artifact pulse arriving before an expected arrival time of a heartbeat pulse; and 

(d) identifying a pulse that arrives at said expected arrival time as a heartbeat pulse only if said artifact pulse 
has subsided. 

21. The method of Claim 19 or 20 wherein step (a) comprises discarding pulses that arrive at a time indicative of a 
change of period between heartbeats greater than approximately 6.75%. 

22. The method of Claim 1 9, 20 or 21 wherein step (b) comprises discarding pulses that have amplitudes varying more 
than 50% from previously received pulses representing heartbeats. 

23. The method of Claim 20 wherein step (c) comprises the steps of: 

after approximately 50% of a period between a last heartbeat pulse and said expected arrival time of a heart- 
beat pulse, setting a noise pulse threshold for detecting artifact pulses; 
lowering said noise pulse threshold until said expected arrival time; and 
identifying non-heartbeat pulses above said noise pulse threshold as artifact pulses. 



PatentansprUche 

1. Vorrichtung zur Bestimmung der Herzrate aus innerhalb eines Korpers erzeugten elektrischen Signalen, umfas- 
send: 

wenigstens zwei elektrische Kontakte (58) zur Messung der elektrischen Signale, wenn sie mit dem Korper 
in Kontakt gebracht werden; 

einen Analogschaltkreis (60), (64), (66), der die elektrischen Signale aufbereitet; 

einen Analog/Digital-Wandler (68), der derart gekoppelt ist, dass er ein analoges Ausgangssignal von dem 
Analogschaltkreis empfftngt und die analogen Signale in eine Vlelzahl von Digitalabtastungen umwandelt; 
einen digitalen Filter zum Empfang der Digitalabtastungen und zurUnterdruckung unerwunschter Signale (70); 
einen Aufbereitungs- bzw. Anreicherungssignalprozessor (92) zum Empfang der gefilterten Daten und zum 
Hervorheben von Signalen in diesem, weiche vorbestimmte Eigenschaften oder QRS-Komplexe In Signalen 
des menschlichen Herzschlags haben, urn angereicherte digitale Daten zu erzeugen; 

dadurch gekennzeichnet, dass der Aufbereitungs- bzw. Anreicherungssignalprozessor aufweist: 

ein Differenzlerglied (96) zur Bestimmung der Steilheit von Spitzen in den gefilterten Daten und zur Erzeugung 
eines Steilheitsslgnals, das die GrdBe und das Vorzelchen der Steilheit eines jeden Berelchs von jeder Spitze 
definiert; 

einen Quadratprozessor (98) zur Quadrierung der Ergebnisse aus dem Differenzierglied; und 
einen Prozessor fur einen beweglichen Durchschnitt (100) zur Berechnung eines beweglichen Durchschnitts 
nur der Signale mit positivem Wert und zum Ausgeben eines Signals des beweglichen Durchschnitts, das den 
beweglichen Durchschnitt Qber die Zeit definiert. 

2. Vorrichtung nach Anspruch 1 , wobel der Analogschaltkreis aufweist: 

einen Differentialverstfirker (60) mit einerVerstarkung von ungefahr 1 - 10 gekoppelt an die elektrischen Kon- . 
takteinhelten zum Verstdrken beliebiger Signale, die von den elektrischen Kontakteinheiten gemessen wer- 
den, und zum Unterdrucken jeglicher allgemeiner Rauschmoden; und 



10 



EP 0 861 045 B1 



einen Analogbandpassfilter (64), der zum Empfang des Ausgangssignals des Differeritialverstarkers ange- 
koppelt ist und ein Durchgangsband von annahernd 5 bis 40 Hz aufweist, 

wobei der Bandpassfilter einen Tiefpass-Analog-Aktivfilter mit einer Abfall-Charakteristlk erster Ordnung mit einer 
5 Eckfrequenz zwischen 25 und 40 Hz und einen Hochpass-Analog-Aktivfilter mit einer Abfall-Charakteristik zweiter 

Ordnung mit einer Eckfrequenz zwischen 5 und 15 Hz aufweist. 

3. Vorrichtung nach Anspruch 1 oder 2, wobei der Analog/Digital-Wandler (68) eine Abtastrate verwendet, die ein 
Vielfaches der AC-Netzspannungsfrequenz ist. 

10 

4. Vorrichtung nach irgendeinem der vorhergehenden Anspruche, wobei der Digitaifilter (90) die Digitalabtastungen 
empfangt und Rauschsignale unterdruckt, die Frequenzen unter5 - 15 Hz haben, und Signale mit Frequenzen 
uber circa 25 - 40 Hz unterdruckt, urn gefilterte Daten zu erzeugen, wobei der Digitalfiiter ein rekursiver Filter mit 
ganzzahligen Koeffizienten ist. 

15 

5. Vorrichtung nach irgendeinem der vorhergehenden Anspruche, wobei das Differenzierglied ein digitaier Signal- 
prozessor ist, der den mathematischen Ausdruck berechnet: 

20 Y^p^ + X^-X^- 2X^)74, 

wobei 

Y n = die Steilheit zu einer beliebigen bestimmten Abtastzeit N darstellt, 
25 2X n = das Doppelte der Amplitude der jungsten Digitaldatenabtastung aus dem Ausgangssignal des vorher- 

gehenden Filter darstellt, 

X n . t = die Amplitude derfoigenden am nachsten gelegenen Digitaldatenabtastung aus dem Ausgangssignal 
des vorhergehenden Filter darstellt, 

X n ^j = die Amplitude der am drittnachsten gelegenen Digitaldatenabtastung in dem Strom von Ausgangssi- 
30 gnalen des vorhergehenden Filter darstellt, und 

2X n -4 = das Doppelte der Amplitude der am viertnachsten gelegenen Digitaldatenabtastung aus dem Aus- 
gangssignal des vorhergehenden Filter darstellt. 

6. Vorrichtung nach irgendeinem der vorhergehenden Anspruche, wobei der Prozessorfur den beweglichen Durch- 
35 schnitt den beweglichen Durchschnitt berechnet als einen Strom von Digitaldatenabtastungen bei dem die am 

n§chsten gelegenen Digitaldatenabtastung den durch Y n dargesteilten beweglichen Durchschnitt darstellt und, 
wobei jeder Y n In dem Datenstrom des beweglichen Durchmessers berechnet wird durch Summation der 8 am 
nfichsten gelegenen Datenproben in dem Strom von Datenproben, der die gefilterten Daten darstellt, und durch 
Division der Summe durch 64. 

40 

7. Vorrichtung nach irgendeinem der vorhergehenden Anspruche, die weiter aufweist einen Signalprozessor zur Be- 
stimmung der Herzrate nach der Verarbeltung, der zur Analyse der Digitaldaten vorgesehen ist, urn daraus die 
Herzrate zu bestimmen. 

4$ 8. Vorrichtung nach Anspruch 7, die weiter aufweist eine Anzeige zum Anzeigen der Herzrate, die durch die Signal- 
prozessoreinheit zur Bestimmung der Herzrate nach der Verarbeitung (109) bestimmt ist wird. 

9. Vorrichtung nach irgendeinem der Anspruche 1 bis 6, die weiter eine Bestimmungseinheit fur den Herzschlag nach 
der Verarbeitung aufweist, die zur Unterscheidung von Herzschlagimpulsen von anderen Impulsen in einer Weise 

so konfiguriert ist, urn; 

Impulse auszusortieren, diezu Zeiten auftreten, die unvereinbarmit einer plausiblen Schlag-zu-Schlag Varia- 
tion der Herzschlagrate 1st, und welche eher Artef akte als HerzschlSge darstellen; und 
Impulse mit Amplituden auszusortieren, die urn mehr als einen vorbestimmten Umfang von vorher empfan- 
55 genen impulsen variieren, welche Herzschlage darstellen. 

10. Vorrichtung nach Anspruch 9, der weiter ausgelegt 1st: 
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zur Identifizierung eines Artefaktimpulses, der vor einer erwarteten Ankunftszeit eines Herzschlagimpulses 
ankommt; und 

zur Identifizierung eines Impulses, derzu der erwarteten Ankunftszeit als ein Herzschlagimpuls ankommt, nur 
falls der Artefaktimpuls abgeklungen ist. 

11. Vorrichtung nach Anspruch 9 oder 10, die zum Aussortieren von Impulsen ausgelegt ist, die zu einer Zeit ankom- 
men, die eine Anderung der Periode zwischen Herzschlagen anzeigt, die grol3er ist als ungefahr 6,75 %. 

12. Vorrichtung nach Anspruch 9, 10 oder 1 1 , die zum Aussortieren von Impulsen ausgelegt ist, die Amplituden haben, 
die urn mehr als 50% von vomer empfangenen Impulsen variieren, welche Herzschlage darstellen. 

13. Vorrichtung nach Anspruch 10, die ausgelegt ist, urn: 

nach ungefahr 50% einer Periode zwischen einem letzten Herzschlagimpuis und der erwarteten Ankunftszeit 
eines Herzschlagimpulses einen Rauschimpulsschwellenwert zur Messung von Artefaktimpulsen zu setzen; 
den Rauschimpulsschwellenwert zu senken bis zu der erwarteten Ankunftszeit; und Nicht-Herzschlagimpulse 
oberhalb des Rauschimpulsschwetlenwerts als Artefaktimpulse zu erkennen. 

14. Verfahren zum Messen einer Herzschlagrate unter unerwunschten Signalen, das die Schritte aufweist: 

Aufnehmen von Signalen von einem Korper durch Verwendung von Elektroden, die mit der Haut des Korpers 
in Kontakt sind; 
VerstSrken der Signale; 

Herausfiltem von Rauschen unterhalb und oberhalb dem Frequenzbereich, in dem die gewiinschte Herzschla- 
grate liegen wird, urn ein gefiltertes Analogsignai vorzusehen; 
Umwandeln des gefilterten Anaiogsignals in eine Vielzahl von Digitalabtastungen; 

digitales Filtern der Digitalabtastungen zum Entfernen unerwunschter Signaie, urn gefiiterte Digitalsignale zu 
erzeugen; 

digitales Verarbeiten der gefilterten Digitalabtastungen zur Aufbereitung bzw. Anreicherung darin enthaltener 
Herzschlagspitzen, urn angereicherte bzw. aufbereitete Digitalabtastungen vorzusehen, gekennzeichnet 
durch 

Differenzierung der gefilterten Digitalabtastungen oben zur Bestimmung der Steilheit der darin codierten Spit- 
zen, urn Dlfferenzierglled-Digitalabtastungen vorzusehen; 

Quadratbildung der Amplitude der in der Differenzierglied-Digitalabtastung codierten Signale, urn quadrierte 
Digitalabtastungen vorzusehen; und 

Berechnung eines beweglichen Durchschnitts aus den quadrierten Digitalabtastungen; und Verarbeitung der 
angereicherten Digitalabtastungen, urn die Herzschlagrate zu bestimmen. 

15. Verfahren nach Anspruch 14, wobei der Schritt des digitalen Filterns der Digitalabtastungen den Schritt aufweist, 
die Probendaten rekursiv zu filtern, die aus den von. dem Korper aufgenommenen Signalen abgeleitet .sind, ent- 
sprechend derfolgenden rekursiven Filterbeschreibung. 

rekursives Tlefpassflltern nach der Glelchung: 

Y n = 1/B(2Y n . 1 -Y lv2 + X n -2X lM + X n . 8 ) 

wobel 

Y n die Ausgangs-Dlgitalabtastung jeder Iteration darstellt, 

und Y^j vorhergehende Ausgangsignale derart darstellt, dass Y n .j das jungste vorhergehende Aus- 
gangssignal darstellt, wahrend Y n . 2 das am zweitnachsten gelegene vorhergehende Ausgangssignal darstellt, 

und wobei 

X Teime wle z. B. X n , X^, X n ^ die vorhergehenden RohdateneingSnge derart darstellen, dass X n den jung- 
sten vorhergehenden Rohdateneingang darstellt, X^ den am viertnachsten gelegenen vorhergehenden Rohda- 
teneingang darstellt und X^q den am achtnachsten gelegenen vorhergehenden Rohdateneingang darstellt, und 

rekursives Bandpassfiltern nach der Glelchung: 



AA 
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Y„ = 1/12 (2V„, • 3Y^ + 2Y^ 3 + X„ - 2X„. e + X fr12 ) 

wobei 

5 

Y n die Ausgangs-Digitalabtastung des Bandpassfiiterungsschrittes bei jeder Iteration darstellt, und 

und Y^g vorhergehende Ausgangs-Digitaiabtastungen des Bandpassfilterungsschrittes derart darstellt, 
dass Y rv1 das jungste vorhergehende Ausgangssignal darstellt, wahrend Y n . 2 das am zweitnachsten gelege- 
nen vorhergehende Ausgarigssignal darstellt, und Y^ das am drittnachsten gelegenen vorhergehende Aus- 
10 gangssignal darstellt, 

und wobei 

XTerme wie z. B. X n , X n ^, X n . 8 die vorhergehenden Probeneingange darstellen, wie etwa Proben abgeleitet 
aus Roh- oder Analogfiiterdaten oder Digitalabtastungsergebnissen einer vorhergehenden Zustandsberechnung, 
is wie dem digitalen Tlefpassfilterungsschritt. 

16. Verfahren nach Anspruch 14, wobei der Schritt der digitalen Filterung weitere Fragmente von Frequenzen uber- 
oder unterhalb des Frequenzbereichs entfernt, in dem der Herzschlag liegt, und.ein Netzgerausch unterdruckt bei 
circa 50 -60 Hz. 

20 

17. Verfahren nach Anspruch 16, wobei der Schritt der Berechnung des beweglichen Durchschnitts die Schritte auf- 
weist: 

Einsatz eines Computers zur Berechnung eines beweglichen Durchschnitts der quadrierten Digitalabtastun- 
25 gen gemaB dem mathematischen Ausdruck, 



wobei, 

35 Y n den beweglichen Durchschnitt zu einer beliebigen bestimmten Abtastzeit N darstellt, 

X N .| die Amplitude der Digitaldatenabtastung zur Zeit N - i darstellt. 

18. Verfahren nach Anspruch 14, wobei die Quadratbildung der Amplitude durch Nachschlagen von Quadratbildungen 
in Nachschlagetabellen durchgefuhrt wird. 

40 

19. Verfahren nach irgendelnem der Ansprtiche 13 bis 18, das zudem das Unterscheiden von Impulse aufweist, die 
durch HerzschlSge von anderen Impulse verursacht sind, durch die Schritte: 

(a) Aussortieren von Impulsen, die zu Zeiten auftreten, die unvereinbar sind mit einer plausiblen Schlag-zu- 
45 Schlag Variation der Herzschlagrate, sondern welche eher Artefakte als Herzschiage darstellen; und 

(b) Aussortieren von Impulsen mit Amplitudes die urn mehr als einen vorbestimmten Umfang von vorher 
empfangenen Impulsen, welche Herzschiage darstellen, variieren. 

20. Verfahren nach Anspruch 19, das zudem die Schritte aufweist: 

so 

(c) Identifizierung eines Artefaktim pulses, der vor einer erwarteten Ankunftszeit eines Herzschlagimpulses 
ankommt; und 

(d) Identifizierung eines Impulses, derzu der erwarteten Ankunftszeit als ein Herzschlagimpuls ankommt, nur 
falls der Artefaktimpuls abgeklungen ist. 

55 

21. Verfahren nach Anspruch 19 oder 20, wobei Schritt (a) das Aussortieren von Impulsen aufweist, die zu einer Zeit 
ankommen, die eine Anderung der Periode zwischen Herzschlagen anzeigt, die grdBer ist als ungefShr 6,75%. 
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22. Verfahren nach Anspruch 19,20 Oder 21 , wobei Schritt (b) das Aussortieren von Impulsen aufweist, die Amplituden 
haben, die urn mehr als 50% von vorher empfangenen Impulsen, welche Herzschlage darstellen, variieren. 

23. Verfahren nach Anspruch 20, wobel Schritt (c) die Schritte aufweist: 

5 

nach ungefahr 50% einer Periode zwischen einem letzten Herzschlagimpuls und der erwarteten Ankunftszeit 
eines Herzschlagimpulses wird ein Rauschimpulsschwellenwert zum Messen eines Artefaktim pulses gesetzt; 
Senken des Rauschimpulsschwellenwert bis zu der erwarteten Ankunftszeit; und 

Identifizierung von Nicht-Herzschlagimpulsen oberhalb des Rauschimpulsschwellenwerts als Artefaktimpuise. 

10 

Revendlcatlons 

1. Dispositif pour determiner !e rythme cardiaque k partir de signaux eiectriques produits k I'interieur d'un corps 
15 comprenant: 

au moins deux contacts eiectriques (58) pour detecter lesdits signaux eiectriques, lorsqu'ils sont places en 
contact avec ledit corps, 

un circuit analogique (60), (64), (66) qui conditionne lesdits signaux electriques; 
20 un convertisseur analogique/numeVique (68) couple de maniere a. recevoir un signal de sortie analogique 

provenant dudit circuit analogique et convertir lesdits signaux analogiques en une plurality d'echantillons nu- 
mdriques; et 

un filtre numerique pour recevoir lesdits echantillons numeriques et supprimer des signaux indesirables (70); 
un processeur de signaux ^intensification (92) pour recevoir lesdites donnees filtrees et mettre en relief, dans 
25 ces donnees, des signaux qui ont des caracteristiques predeterminers ou des complexes QRS predetermines 

dans des signaux de battements cardiaques humains de maniere k produire des donnees num6riques ampli- 
fies, 

caracterise en ce que ledit processeur de signaux ^intensification comprend: 

30 

un circuit differenciateur (96) pour determiner la pente de pics dans lesdites donnees filtrees et produire un 
signal de pente qui definit ^amplitude et le signe des pentes de chaque partie de chacun desdits pics; 
un processeur d'6l6vation au carre (98) pour 6lever au carre des resultats fournis par ledit circuit differenciateur; 
et 

35 un processeur (100) de formation de moyenne mobile pour calculer une moyenne mobile desdites valeurs 

positives uniquement du signal et deiivrer un signal de moyenne module qui definit ladite moyenne mobile 
dans !e temps. 

2. Dispositif selon la revendication 1 , dans iequel ledit circuit analogique comprend: 

40 

un ampiificateur differential (60) possedant un gain compris entre environ 1 et 10 et couple auxdits moyens 
de contact eiectrique pour amplifier tous les signaux detectes par lesdits moyens de contact eiectrique et 
supprimer tout bruit en mode commun, et 

un filtre passe-bande analogique (64) couple de maniere k recevoir le signal de sortie dudit ampiificateur 
45 differentiel et possedant une bande passante comprise entre environ 5 et 40 hertz, ledit filtre passe-bande 

comprenant un filtre passe-bas analogique actif possedant une caracteristique de coupure du premier ordre 
avec une frequence limite entre 25 et 40 hertz, et un filtre passe-haut analogique actif possedant une seconde 
caracteristique de coupure avec une frequence limite comprise entre 5 et 15 hertz. 

so 3. Dispositif selon la revendication 1 ou la revendication 2, dans Iequel ledit convertisseur analogique/numerique 
(68) utilise une frequence d'echantillonnage qui est un multiple de la frequence de la tension alternative de la ligne 
d'alimentation. 

4. Dispositif selon Tune quelconque des revendlcatlons prec6dentes, dans Iequel ledit filtre numerique (90) regoit 
55 lesdits echantillons num6riques et supprime des signaux de bruit ayant des frequences Inf6rieures k 5-15 hertz 
et des signaux possedant des frequences superieures k environ 25-40 hertz pour produire des donnees filtrees, 
ledit filtre num6rique etant un filtre recurslf possedant des coefficients formes de nombres entiers. 
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5. Dispositif selon Tune quelconque des revendications pr6cedentes ( dans lequel ledit circuit difterenciateur est un 
processeur de signaux num6riques qui calcule ('expression math6matique 

5 Y^iaX^X^-X^- 2X^4, 

dans laquetle 

Y n = repr6sente la p.ente k n'importe quel instant particulier d'echantillonnage N, 
10 2X n = represente le double de I'amplitude de I'echantillon de donnees numeriques le plus recent fourni par la 

sortie du filtre precedent, 

X n>1 = represente Pamplitude de i'echantillon de donnees num6riques le plus recent en second, fourni par la 
sortie du filtre precedent, 

X n . 3 = represente I'amplitude du troisieme 6chantillon de donnees num6riques le plus recent dans ie flux 
is deiivr6 par la sortie du filtre precedent, et 

2X n . 4 = represente le double de ['amplitude du quatrfeme echantillon de donnees numeriques le plus recent 
deiivre par la sortie du filtre precedent 

6. Dispositif selon Tune quelconque des revendications precedentes, dans lequel ledit processeur de formation de 
20 la moyenne mobile calcule ladite moyenne mobile sous la forme d'un flux d'6chantillons de donnees numeriques, 

I'echantillon de donnees numeriques le plus recent representant la moyenne mobile representee par Y N , et dans 
lequel chaque Y N dans ledit flux de donnees de la moyenne numerique est calcule par sommation des 8 echan- 
tillons de donnees les plus r6cents dans ledit flux d'6chantillons de donn6es representant lesdites donn6es filtr6es 
et division de la somme par 64. 

25 

7. Dispositif selon I'une quelconque des revendications precedentes, comprenant en outre un processus de signaux 
de determination du rythme cardiaque aprfcs traitement, agence pour analyser lesdites donnees numeriques pour 
determiner, k partir de \k, le rythme cardiaque. 

so 8. Dispositif selon la revendication 7, comprenant en outre un dispositif d'affichage pour aff icher ledit rythme cardiaque 
determine par les moyens (109) formant processeur de signaux de determination du rythme cardiaque aprfcs 
traitement. 

9. Dispositif selon I'une quelconque des revendications 1 k 6, comprenant en outre des moyens de determination du 
35 rythme cardiaque aprfcs traitement configures de manifcre k etablir une distinction entre des impulsions de batte- 

ments cardiaques vis-6-vis d'autres impulsions, par le fait que ces moyens sont agences pour: 

eiiminer des impulsions qui arrivent k des instants incompatibles avec une variation plausible d'un battement 
au suivant dans le rythme cardiaque, et ce comme repr6sentant des artefacts plutdt que des battements 
to cardiaques; et 

eiiminer des impulsions ayant des amplitudes variant de plus d'une vateur predeterminee par rapport k des 
impulsions recues precedemment et representant des battements cardiaques. 

10. Dispositif selon la revendication 9, configure en outre de manidre k: 

45 

identifier une impulsion d'artefact arrivant avant un Instant d'arriv6e attendu d'une impulsion de battement 
cardiaque; et 

identifier une impulsion qui arrive audit instant d'arrivee attendu en tant qu'impulsion de battement cardiaque 
uniquement si ladite impulsion d'artefact a disparu. 

so 

11. Dispositif selon la revendication 9 ou 10, agence de manifcre k eiiminer les impulsions qui arrivent k un instant 
indicatif d'un changement de p6riode entre des battements cardiaques, superieur k environ 6,75 %. 

1 2. Dispositif selon la revendication 9, 1 0 ou 1 1 , agence de manifcre k eiiminer des impulsions qui ont des amplitudes 
55 variant de plus de 50 % par rapport k des impulsions recues anterieurement et representant des battements 

cardiaques. 

13. Dispositif selon la revendication 10, agence de manifere k\ 
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au bout d'environ 50 % d'une pEriode entre une derniEre impulsion de battement cardiaque et ledit instant 
d'arrivEe attendu d'une impulsion de battement cardiaque, regler un seuil d'impulsion de bruit pour detecter 
des impulsions d'artefact; 

reduire ledit seuii d'impulsion de bruit jusqu'audit instant d'arrivee attendu; et 
5 identifier des impulsions qui ne sont pas des battements cardiaques, au-dessus dudit seui! d'impulsion de 

bruit, en tant qu'impulsions d'artefact. 

14. Procede pour detecter un rythme cardiaque parmi des signaux indesirables, comprenant les etapes consistant 

10 detecter des signaux delivres par un corps en utilisant des electrodes placees en contact avec la peau du corps; 

amplifier les signaux; 

eliminer par filtrage un bruit situe au-dessous et au-dessus de la gamme des frequences, dans laquelle le 
rythme cardiaque se situe, de maniere & foumir un signal analogique filtre; 
convertir le signal analogique filtre en une plurality d'echantilions numeriques; 
15 fiitrer numeriquement les Echantillons numeriques pour eliminer des signaux indesirables afin de produire des 

signaux numeriques filtres; 

traitor numeriquement les echantillons numeriques filtrEs pour intensifier des pics de battements cardiaques 
represented dans ces echantillons de maniere k fournir des echantillons numeriques amplifies, 

20 caracterise en ce 

qu'on differentie les Echantillons numeriques filtres pour determiner la pente de pics codes dans ces echan- 
tillons pour fournir des echantillons numeriques du circuit differenciateur, 

qu'on Eleve au carre I'amplitude de signaux codes dans ledit Echantillon numErique du circuit differenciateur 
25 pour produire des echantillons numeriques ElevEs au carrE, et 

qu'on calcule une moyenne mobile desdits Echantillons numeriques eleves au carre, et qu'on traite lesdits 
echantillons numeriques amplifies pour determiner le rythme cardiaque. 

15. Procede selon la revendication 14, selon lequel I'etape de filtrage numErique des Echantillons numeriques com- 
30 prend Petape de filtrage rEcursif des donnEes Echantillons dErivEes des signaux dEtectEs & partir dudit corps con- 

formEment aux spEcifications suivantes de filtre rEcursif , avec filtrage passe-bas rEcursif conformEment a. ['Equa- 
tion: 

35 Y n = 1/8(2Y n . l -Y rv2 + X n -2X Iv4 + X n . 8 ) 

dans laquelle 

Y n reprEsente I'Echantillon numErique de sortie lors de chaque iteration, 
40 Y n . 1 et Y n . 2 reprEsentent diffErents signaux de sortie, par exemple que Y n . t reprEsente le signal de sortie 

precedent le plus recent, tandis que Y^ reprEsente le signal de sortie precedent le plus recent en second, 

et selon lequel 

les termes X tels que X n , X n ^ et X^ reprEsentent des entrEes de donnEes brutes rEcentes, par exemple X n 
45 reprEsente PentrEe de donnEes brutes la plus recente, X^ reprEsente I'entrEe de donnEes brutes la plus rEcente 

en 4-Eme position et X n . 8 reprEsente I'entrEe de donnEes brutes la plus rEcente en 8-eme position, et 
et filtrage passe-bande rEcursif conformEment k I'Equation: 

50 Y n - 1/iaPYn-l * 3Y n-2 + 2Y n . 3 + X n - 2X^3 + X n _ 12 ) 

dans laquelle 

Y n reprEsente I'Echantillon numErique de sortie fourni par I'Etage de filtrage passe-bande lors de chaque itE- 
55 ration, et 

Y n-1 et Y n . 2 reprEsentent des echantillons numeriques de sortie precedents de I'etape de filtrage passe-bande 
de telle sorte que Y rv1 reprEsente le signal de sortie precedent le plus recent, tandis que Y rv2 reprEsente le 
signal de sortie precedent le plus rEcent en second, et Y n . 3 reprEsente le signal de sortie prEcEdent le plus 
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recent en troisieme position, 
et selon lequel 

les termes X tels que X n> Xp. 4 et X^ repr6sentent des entrees d'6chantillons recentes tels que des echantillons 
5 tires des donnees filtr£es brutes ou analogiques ou des resultats d'echantillons numeriques du calcul d'un etage 

precedent, tel que l'6tape de filtrage passe-bas numerique. 

16. Proced6 selon la revendication 14, selon lequel ladite etape de filtrage numerique eiimine d'autres restes de fre- 
quences au-dessus et au-dessous de la gamme de frequences, dans laquelle ledlt rythme cardiaque se situe, et 

10 sert a supprimer le ronflement dans la llgne d'alimentation a environ 50-60 hertz. 

17. Procede selon la revendication 16, selon lequel ladite etape de calcul de ladite moyenne mobile comprend les 
Stapes consistent a: 

15 utiliser un ordinateur pour calculer une moyenne mobile desdits echantillons numeriques eiev6s au carre 

conformement a Pexpression mathematique: 

1 7 
64 i«0 

25 dans laquelle 

Y N represente la moyenne mobile pour n'importe quel instant particulier d'echantillonnage N, 
x n-i repr6sente I'amplitude de rechantillon de donnees numeriques a I'instant N-i. 

30 18. Procede selon la revendication 14, selon lequel ladite elevation au carre de I'amplitude est executee par consul- 
tation de carres dans une table de consultation. 

19. Proc6d6 selon Tune quelconque des revendications 13 a 18, comprenant en outre une operation servant a distin- 
guerdes impulsions provoquees paries battementscardiaques par rapport a d'autres impulsions, moyennant les 

35 etapes consistant a: 

(a) eiiminer des impulsions qui arrivent a des instants incompatibles avec une variation plausible d'un batte- 
ment au suivant dans le rythme cardiaque, et ce comme repr6sentant des artefacts plut6t que des battements 
cardiaques; et 

40 (b) eiiminer des Impulsions ayant des amplitudes variant de plus d'une valeur predeterminee par rapport a 

des impulsions regues precedemment et representant des battements cardiaques. 

20. Procede selon la revendication 1 9, mettant en oeuvre les etapes consistant a: 

^5 (c) identifier une impulsion d'artefact arrivant avant un instant d'arrivee attendu d'une impulsion de battement 

cardiaque; et 

(d) identifier une impulsion qui arrive audit instant d'arrivee attendu en tant qu'impulsion de battement cardia- 
que uniquement si ladite impulsion d'artefact a disparu. 

50 21. Proc6de selon la revendication 19 ou 20, selon lequel retape (a) comprend reiimination d'impulsions qui arrivent 
a un instant indicatif d'une variation de la p6riode entre des battements cardiaques, superieure a environ 6,75 %. 

22. Proc6d6 selon la revendication 19, 20 ou 21, selon lequel retape (b) comprend reiimination d'impulsions qui ont 
des amplitudes variant de plus de 50 % par rapport a des impulsions revues pr6c6demment et representant des 

55 battements cardiaques. 

23. Proc6d6 selon la revendication 2, selon retape (c) comprend les etapes a: 
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au bout d'environ 50 % d'une p6riode entre une dernifcre impulsion de battement cardiaque et ledit instant 
d'arriv6e attendu d'une impulsion de battement cardiaque, r6gler un seuil d'impulsion de bruit pour d6tecter 
des impulsions d'art6fact; 

r6dulre ledit seuil d'impulsion de bruit jusqu'audit instant d'arriv6e attendu; et 
s identifier des impulsions qui ne sont pas des battements cardiaques, au-dessus dudit seuil d'impulsion de 

bruit, en tant qu'impulsions d'artefact. 
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